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ABSTRACT 


The  following  model  problem  is  studied; 


ft 


r  1L  fr.g  .  JL  <  A  iHn 
r  3r  (r3  3r  3z  (8  3z)] 


f 


ri  :  u  «  0 


where  ft  is  a  bounded  open  domain  with  r  <  0  in  (r,z)  plane,  I*  *  3ft\To, 
To  *  3ft  n  {(r,z)  :  r  *  0}  •  We  introduce  weighted  Sobolev  spaces  V*(k  * 
1,2),  and  prove : 

1  2 

(1)  The  problem  has  a  unique  solution  u,  and  u  6  VQ(ft)  n  V  (ft). 

(2)  The  linear  finite  element  solution  u^  exists  and  is  unique. 

(3)  The  error  u-uh  in  "energy  norm"  is  of  0(h).  Particularly,  if  ft 

is  a  polygon,  then 

lu  -  Vi,a  '  olh) 

-  Vo, a  '  0(h2> 

k 

where  1**^  ^(k  *  1,2)  are  the  V  norms. 
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Introduction 

The  numerical  solution  of  singular  boundary  value  pobblems/fiave  been  studied  by 

/  / 

several  authors.  The  finite  difference  methods  and  its  theory  for  a  type  of  two- 
dimensional  singular  boundary  value  problems  are  given  in  [ID],  [13].  The  finite  element 
method  for  axisymmetrlc  elastic  solid  is  proposed  in  1)6}.  15},  [11),  [14}  and  [20),  gives 

a  proof  of  the  convergence  of  the  finite  element  methods  for  one  dimensional  singular 
problems.  [12]  proves  the  "optimal"  order  of  convergence  jf or  the  method  of  [16)  provided 
the  loads  are  axlsymmetric  and  the  solution  is  in  C  (ft).  The  convergence  of  the  linear 
finite  element  method  for  two  dimensional  singular  Dirichlet  problem  is  proved  in  [18).  In 
this  paper  we  will  prove  the  so-called  "optimal"  order  of  convergence  of  the  linear  firu^e 
element  method  for  the  following  model  problem: 


ft  : 


r  1  3  r  a  3u  i  3  f  *  3u  1 1 


r1!  u 


r  3r 

0 


(1.1) 


where  ft  is  a  bounded  open  domain  with  r  >  0  in  (r,z) -plane, 

T1  *  3ft/ rQ,  «  3ft  n  {(r,z):  r  »  o}. 

We  assume: 

(i)  The  function  3  is  uniformly  Lipschitz  continuous  in  ft. 

(ii)  0  >  6q  >  0,  3q  is  a  constant. 

(iii)  t/2i  e  l2{0). 

f' 
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Under  the  (x,y,z)  coordinate  system  we  have 


,  •  9**  -i  ,  Jvi  3vi  i 

-  J.  vn  57  r  dxdyd*  -  J,  (g—  CO.  e  +  5 —  .in  0)  **r  dxdydr 

0  a  y 


(2.5) 


— \  . ]  da*  *  *  i 

Since  r  **  and  r  ^  are  bounded  in  U  ,  v  e  H  (8*),  we  may  take  the  limit  through 

(2.5)  as  n  ♦  •  and  hence  we  obtain  (2.5)  ae  well  as  (2.4)  with  v,v*  replacing  vR,  vR* 


respectively.  (2.1)  is  proved. 

Simple  calculation  derives  the  following  results. 


Q.E.D. 


Corollary  2.1.  If  v*  e  H  (8*),  then 


3v#  3v 

K  *  Ti 


cos  6 


3y  sin© 
36  — 


iv*  3v  A  3y  cos© 

~  “  8in  0  +  ao 


in  8* 


3y  3r  d©  r 

Corollary  2.2.  Assume  v  is  independent  of  0.  Then  we  have  for  v*  e  H^Q*): 


for  v*  e  H  (8*) 


3v*  3v  Q  3v*  3v  .  a 

j;  *  57 cos  6 '  a?  ■  a?  *ln  9 » 


a2v* 

a*2 


32v  2.  3v  sin29 

— r  cos  0  +  *5 - 

3r2  3r  r 


.in20  + 


3v  cos20 
3r  r 


32y* 

3x3y 


1 

r  3r} 


sin  0  cos  0  , 


32v*  32v  32v*  32v  32v*  32v 

3z2  3z2  *  ^r3z  cos  9  dydz  3r3z 
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3.  Spaces  V1.  V2([2j,  [17] ) 

We  define  functionals  l*loK“0#1#2#  as  follows: 


lvln  o  “  <  /  v2rdrdzr2 


,V,10‘(  I  ,»<‘v,o  a>^ 

’*  l«l<l  °' 


Ivl,  Q  -  (  l  I9°vl^  Q  +  I  1  I2  l„  -)1/2 
2.0  |oj<2  0,0  r  3r  0,0 


Definition  3.1.  Assume  that  D  is  an  open  or  closed  set  in  (r ,z) -plane,  D*  the 
correspondent  axisymmetric  set  in  (x,y,z) -space,  A(D>  the  set  of  real  functions  defined  in 
D,  and 

A*(D*)  -  {v*  :  v*  real  function  defined  in  D*,  and  there 

exists  v  e  A(D)  such  that  v*(x,y,z)  »  v ( J x‘ 2  +  y2 , z ) } • 

We  define  a  mapping  T:  A*(D*)  ♦  A(D)  as  follows: 

Tv* (x,y, z)  -  v( r,z) 

Obviously,  the  mapping  T  is  one-to-one* 

Definition  3.2.  Uk(0*)  -  H  fo*)  nA»(0»),  k  -  0,1,2. 

k  lr 

It  is  easy  to  see  that  U  (ft*)  is  a  closed  subspace  in  H  (ft*).  Now  establish  the 

relations  between  the  norms  1*1  and  the  functionals  I -Ik, ft  for  the  elements  of 

k  n  (ft* ) 

U  (ft*). 


3. 1 .  Assume  v*  €  U  (ft*),  v  -  Tv*.  then  Ivl^  ^  <  •,  and 


lv*l 


*.2 


.  -  2»lvl2  Q  V  u*  e  U  k(ft*),  k  -  0,1  . 

H  (ft* ) 


(3.1) 


T“  lvl9  Q  <  lv#|2o  K  2llvl2  ¥  u*  e  U  2(ft*) 
2  2'M  H2 ( ft* )  2'"' 


(3.2) 


Proof.  By  direct  computation  and  corollary  2.2. 
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Definition  3.3.  Vk(fl)  -  {vjv  -  Tv*,  v*  e  (J  k(fl*)},  k  -  0,  1,  2. 

It  follows  froa  leans  3.1  end  the  closeness  of  Uk(fl*)  in  Hk(ft*)  that 

V  (ft),  k  -  0,  1,  2,  are  Banach  apace*.  We  need  the  following  aubapace 
v’<0)  of  v1 (0) * 

v’(fl)  -  {viv  -  TV*,  V*  c  u  Vo*)  n  h’(0**}. 

Let  v  9  vj(fl),  v  "  Tv*,  tr  v*  be  the  trace  of  v*  on  30*.  He  define  T(trv*)  as 
the  trace  of  v  on  P^,  Obvioualy,  it  la  aero. 

By  lrnna  3.1  and  the  embedding  theorem*  of  H  (ft*).  We  obtain  the  correspondent 
theorem*  of  V*(ft*).  Particularly,  we  have  the  following  result. 

Lemma  3.2.  There  exist*  a  constant  C1  such  that 


Finally ,  the  following  statement  on  denseness  may  be  proved  (see  [17]  for  V1 .  the 
proof  is  similar  for  V2). 

Lemma  3.3.  Assume  that  the  domain  ft  has  a  locally  Lipschita  Boundary.  Then  C  (ft) 
is  dens,  in  vNfl),  k  -  1,  2. 

Remark  3.1.  Lemma  3.3  is  not  a  direct  corollary  of  the  denseness  theorem  of 

k  _  tr  •  — 

H  (ft*).  If  v*  C  H  (ft*),  then  there  exists  a  sequence  v^*  e  C  (ft*)  converging  to 

v*  in  hNq*).  But  we  can  not  claim  that  v  *  e  A*(ft*). 

n 

Remark  3.2.  The  facts  v  e  C  (ft)  and  v  «  Tv*  do  not  imply  that  v*  e  C**(ft*). 

Counter  example:  v  -  r.  But  v  e  C°(ft)  <“>  v*  e  C°(ft*). 


4.  Solution  of  problem  (1*1) 

We  define  a  bilinear  form  B(*,*)  on  V*  (SI)  x  V^Q)  and  a  linear  functional 
F(* )  on  V^(Q)  as  follows; 

B(u,v)  -  /  B(-|“  I7  +  1^  I^)rdrd* 

F(v)  “  /  fvrdrds 
Q 


Then  we  have  the  variational  formulation  of  problem  (1.1  )x  Find  u  e  vJ(Q)  such  that 

B(u,v)  -  F(v),  V  V  e  v’(0)  (4.1) 

From  now  on  we  assume  that  U  has  a  locally  Lipschitz  boundary. 

Theorem  4 . 1  Problem  (4.1)  has  a  unique  solution. 

Proof;  It  follows  from  lemma  3.2  and  assumptions  (i)-(ii)  that  the  bilinear  form  B(u,v) 
is  coercive  and  continuous  on  V^(fl)  x  vJ(Q).  And  the  linear  functional  F(v)  is 
continuous  on  V^(Q)  by  virtue  of  assumption  (iii).  Hence  the  conclusion  of  the  theorem 
is  a  result  of  the  Lax-Milgram  theorem.  Q.E.D. 

Remark  4.1.  Let  u  be  the  solution  of  problem  (4.1).  Since  B(u.v)  is  symmetric, 
U  is  also  the  solution  of  the  following  problem;  Find  u  €  vjJ(Q)  such  that 

J(u)  -  min  J(v) 

v  e  v’(Q) 


where  J(v)  *  B(v,v)  -  2F(v). 

Consider  the  boundary  value  problem  in  Q*  corresponding  to  problem  (1.1); 


dfl#:  w#  *  0 


(4.2) 


where  0#  *  T  ^3,  f*  -  T  1f.  Correspondent  variational  problem  is;  Find  w#  e 
such  that 

«  r^v*),  ¥  v*  e  h’(0*)  (4.3) 
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where 


a  dw*dv 

V"**v*>  -  JZ  *  ^  ^  +  T,  37  )dxdyd' 


F.tv*)  -  /  f*v*dxdyds 
1  ft* 


Prom  now  on  we  assumes 

(iv).  The  boundary  3ft*  is  smooth  snough  to  ensura  that  problem  (4*3)  has  a  unique 
1  2 

solution  w*  and  w*  *  Hg(ft*)  n  H  (ft*)*  For  example,  we  may  assume  that  3ft*  is  of 
class  C2  (sea,  for  instance,  [9,  p.176])  or  that  the  domain  ft*  is  convex. 

Theorem  4.2.  Let  u  be  the  solution  of  problem  (4.1).  Then 


u  e  v*(ft)  n  v2(ft) 

Proof i  Let  u*  •  T-1u.  We  define  for  v*  e  Hg(ft*>  that 

v(r,9,s)  -  v*(x,y,s) 

_  2W 

v(r,x)  -  /  v(r,6,s)d0 

0 


It  is  easily  proved  that 

v  €  V*(ft) 


(4.4) 


3v 

3r 


2*. 

1  £*>' 


3v 

3z 


2wa 

i 

o 


(4.5) 


Now  we  prove  that  u*  is  the  solution  of  problem  (4.3).  It  follows  from  lemma  2.1,  (4.4) 
(4.5)  and  (4.1)  that  for  v*  e  H ’(»*) 

BjtuSv*)  -  P^v*)  -  -  fu)d9]rdrdi 

ft  o 

-  /  ( <§7  §7  ♦  f^>  -  fv]rdrd2  -  B(u,v)  -  F(v)  -  0. 

H.nc.  u*  1.  the  solution  of  problan  (4.3),  and  u»  e  h’(Q*)  H2(fl*)  by  assumption 

(iv).  According  to  definition  3.3  we  obtain  the  conclusion  of  the  theorem.  Q.E.D. 


5.  Linear  finite  element  solution  and  its  order  of  convergence  order 

Assume  that  the  domain  Q  is  convex.  Let  T.  *  {c<l,...,Cm}  be  a  triangula 

h  I 

tion  of  Q,  h^  the  maximum  edge  of  the  triangle  C^,  ©^  the  minimum  angle  of 

m 

C  ,  h  -  max  h  ,  6  -  min  0  ,  Q  -  l  C  ,  assume: 

1  i  1  1  n  i-1 

(a).  ©  >  ©  >  0,  ©  is  independent  of  h  ([7],  [19]). 

o  o 

Define  a  linear  finite  element  space  Vh  as  follows: 


v*1  -  { v.  C  C^(Q)  :  v.  is  linear  function  in  C.  ,  i  ■ 
h  n  i 

vh  -  0  on  (Q  -  Qh)  u  Tl}  # 

Then  it  is  easy  to  prove  that  V11  C  vJ(Q).  We  have  the  correspondent  discrete  problem  for 

problem  (4.1):  Find  u^  8  V*1  such  that 

B(uh#vh)  -  F(vh),  V  vh  e  Vh  (5.1) 

Theorem  5.1.  Problem  (3.1)  has  a  unique  solution. 

The  proof  is  similar  to  that  of  theorem  4.1. 

Remark  5.1.  The  solution  Uh  of  (5.1)  is  also  the  solution  of  the  minimiration  problem: 

Find  li  e  Vh  such  that  J(U.  )  -  min  J(V.  ). 
h  h  vh  e  ^  h 

Assume  that  u  is  the  solution  of  problem  (4.1),  Uj  the  piecewise  linear 

intrpolation  corresponding  to  the  triangulation  T^.  For  any  triangle  C  e  T^,  we  now 

estimate  lu  -  u  I 4  .  Let  Pj  -  (rH,rH),  j  -  1,2,3  be  the  vertexes  of  C, 

I  1 ,c  J  J 

A  (r,s),j  ■  1,2,3  the  so-called  barycentric  coordinates  ([4,  p.  45]),  i.e.  the  basis 
functions  for  the  linear  interpolation  on  C: 


Then  we  have  for  any  function 


Particularly, 


vv 


ij 


v  defined  on 


l  X  (P)V(P.) 

J  3  3 


1.  I  Xjfj  -  r. 


(i,j  -  1,2,3) 

c  and  its  linear  interpolation  V^: 

«  vJ(p),  v  p  e  c, 

A  *  -  r,  V(  r, z)  e  C, 

j  3  3 


(5.2) 


(5.3) 


It  follows  from  (5.3)  that 

l  A  (r  -  r)  -  l  X  (z  -  *)  -  0,  ¥  (r,z)  e  C  (5.4) 

j  3  3  j  3  3 
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Integrating  by  parts  the  integrals  in  the  second  sum,  noting  (5.4)  and  that 


at  “V'V1 


djV(Mj  ^ ' 


we  derive  from  (5.8)  that 


3vt  a  1  ax,  2 
»T-  37’  \{  (1  -  ‘>3^  V(Mj>dt 


It  follows  from  the  uniform  basis  condition  (a)  that  (see,  for  instance,  [8]  or  [15, 


P* 137]  ) 


<  l^h 


(5.10) 


where  h  is  the  maximum  edge  of  c,  M  *  4/sin  6  •  Hence  we  have 

l  o 


—  I 
ar1 


<  M^h 


l  /  (1-t)|d*v(M  )|dt  , 
j  0  33 


and  then 


0V  .  . 

/  lj-r  -  f^|2rdrdz  <  M2h"2  /  (  \  f\l  -  t/4  (1  -  t)5/4|d2v(M  )|dt)2rdrdz 
c  r  cjo  33 

<  3M2h“2  l  j  (/’(I  -  t)5/2|d2v(M.)|2dt*/1(1  -  t)“  1/2  dt)rdrdz 
j  c  0  33  0 


=  6M2h-2  l  J1 dt  /  (1  -  t)5/2|[(r  -  r  +  U  -  z)4— ]  2v(M  . )  |  2rdrdz 

1  j  0  c  3  r  3  *  3 


32v(M. ) 


32v(M , ) 


32v(M . ) 


0  5  1  C  /  J  A  X  A  '  '  -i  '  . 

<  6M2h_2E  /  at  /  (1  -  t)5/2h4(| - 1|  +  2|gr3z  J  I  +  | — Y-1  l>  rdrd* 

j  0  c  3r3z  3z 


2 

where  M 2  *  72  M  ^  •  Hake  variable  transformations  in  the  integrals  as  follows: 

5  “  V*  r(1-t)  ,  n  -  Zj  ♦  z(1-t) 

Then  *  (C#n),  and  the  triangle  C  reduces  to  a  similar  triangle  C^,t  with  the 
similarity  transformation  center  P^«  Hence  the  right  side  of  (5*11)  becomes: 


Since 


we  have 


%  -  vh  -  o  on  a  -  ah  , 


Vu  •  Vu  •  V  “  1".h  Vu  -  vh'u  "  V 


v.  e  v 
n 


Define  u  ■  0  on  8  -  0  .  Then  u  e  V 
i  n  I 


So 


Vu  •  Vu  •  V <  8h(u  •  UI'U  ■  V  *  “x  ®',u  - 


(5.14) 


By  virtue  of  lenma  5. 1  we  have 


a  m 

,u  -  UT*1  o  “I  ,u  -  ut'i  *  Mh2  l  ,u,25  *  Mh2lul2,  o  (5.15) 

1  1#0h  i-1  1  1/Ci  i-1  2,Ci  2#Q 


(5.14)  and  (5.15)  prove  that  (5.13)  is  valid.  Q.E.D. 

*  8,  B^CVfV)  *  B(v,v).  Since  B(u,v)  is  coercive  on 

VjJ(fl),  we  have 

Corollary  5.1.  If  8  is  a  polygon,  then 

,u  •  Vi, a  "  0(h> 

,u  •  Vo,a  “  0(h2) 
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